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^ • Charge transport in the diffusive normal metal (DN) / insulator / s- and d-wave superconductor 

h^ ' junctions is studied in the presence of magnetic impurities in DN in the framework of the quasi- 

^H , classical Usadel equations with the generalized boundary conditions. The cases of s- and d-wave 

(— ^ ■ superconducting electrodes are considered. The junction conductance is calculated as a function of 

fT^ ' a bias voltage for various parameters of the DN metal: resistivity, Thouless energy, the magnetic 

impurity scattering rate and the transparency of the insulating barrier between DN and a super- 
conductor. It is shown that the proximity effect is suppressed by magnetic impurity scattering in 
DN for any value of the barrier transparency. In low-transparent s-wave junctions this leads to 
the suppression of the normalized zero-bias conductance. In contrast to that, in high transparent 
junctions zero-bias conductance is enhanced by magnetic impurity scattering. The physical origin 
5— ( ' of this effect is discussed. For the d-wave junctions, the dependence on the misorientation angle 

t^' a between the interface normal and the crystal axis of a superconductor is studied. The zero-bias 

^ , conductance peak is suppressed by the magnetic impurity scattering only for low transparent junc- 

tions with a ~ 0. In other cases the conductance of the d-wave junctions does not depend on 
"t2 ' the magnetic impurity scattering due to strong suppression of the proximity effect by the midgap 

Andreev resonant states. 
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I. INTRODUCTION 



o, 

CJ , Nowadays, thanks to the nanofabrication technique, detailed experimental studies of the electron coherence in 

mesoscopic superconducting systems become possible, where the Andreev reflection^'^''^ plays an important role in 

p\| ', the low energy transport. In diffusive normal metal / superconductor (DN/S) junctions, the phase coherence between 

I> • incoming electrons and Andreev reflected holes persists in DN at a mesoscopic length scale and results in strong 

l^ interference effects on the probability of Andreev reflectiorki. 

[^ : One of the remarkable experimental manifestations of the coherent Andreev reflection is the zero bias conductance 

peak (ZBCP) in DN/S junctions^'Q''^-^-^-^°-"'^^i^^i^'^i^^. The physics of ZBCP was extensively studied theoretically using 

scattering matrix approach -^^i^'^i-^^i^^i^°i^-^ and the quasiclassical Green's function techmquo^^i^^i^^i^'''i^^i^^i^°i'^-^i'^^i^^i^^. 

■^ Volkov, Zaitsev and Klapwijk (VZK)— explained the origin of the ZBCP in DN/S junctions in the framework of the 

C^ quasiclassical theory by solving the Usadel equations^'^ with the Kupriyanov and Lukichev (KL) boundary condition 

^ for the Keldysh-Nambu Green's functionSi. According to the VZK theory the ZBCP is due to the enhancement of the 

Cd _ pair amplitude in DN by the proximity effect. The influence of the magnetic impurity scattering on the bias voltage 

^ ' dependent conductance was also studied within this approacb22i2L2^. 

I ■ Recently the VZK theory for s-wave superconductors was extended by Tanaka et ali^ using more general boundary 

conditions provided by the circuit theory of Nazaro"\i2&. These boundary conditions treat an interface as an arbitrary 
connector between diffusive metals. The connector is characterized by a set of transmission coefficients ranging from 
^ , a ballistic point contact to a tunnel junction. The boundary conditions coincide with the KL conditions when a 
connector is diffusive or transmission coefhcients are low, while the BTK theory^ is reproduced in the ballistic regime. 
The extended VZK theorji^Lii revealed a number of new features like a U -shaped gap like structure and a crossover 
from a zero bias conductance peak (ZBCP) to a zero bias conductance dip (ZBCD). These phenomena are relevant for 
5— i the actual junctions in which the barrier transparency is not necessarily small. However, the influence of the magnetic 

impurity scattering in DN on the charge transport was not studied in this regime. 

The generalized VZK theory was recently applied also to unconventional superconducting junctions i^iia. The 
formation of the midgap Andreev resonant states (MARS) at the interface of unconventional supcrconductor822i2242iii 
is naturally taken into account in this approach^^*^. It was demonstrated that the formation of MARS in DN/d- 
wave superconductor(DN/(i) junctions strongly competes with the proximity effect. Remarkable recent advances in 
experiments on tunneling in high Tc cuprates'^^ stimulate an interest to the problem of an influence of the magnetic 
impurity scattering on a charge transport in DN/d junctions. 

In the present paper the generalized VZK theory is applied to the study of an influence of the magnetic impurity 
scattering in the DN on the conductance in DN/S where S is either s- or d-wave superconductor. The parameters 



of the problem are the height of the insulating barrier at the DN/S interface, the resistance Rd, the magnetic 
impurity scattering rate 7, the Thouless energy Exh in DN and the angle a between the normal to the interface 
and the crystal axis of d-wave superconductors. We shall focus on the dependence of the normalized conductance 
axieV) — as{eV)/aNieV), on the bias voltage V, where as^]^\{eV) are the conductances in the superconducting 
(normal ) state. The organization of the paper is as follows. In section II the detailed derivation of the expression 
for the normalized conductance is provided. In sections III the results of calculations of axieV) are presented for s- 
and d-wave junctions separately and physical explanation of the results is given. In section IV the summary of the 
obtained results and the conclusions are presented. In this paper we restrict ourselves to the low-temperature regime 
T « Tc and put ks — ft — 1. 

II. FORMULATION 

In this section we introduce the model and the formalism. We consider a junction consisting of normal and 
superconducting reservoirs connected by a quasi-one-dimensional diffusive conductor (DN) with a length L much 
larger than the mean free path. This structure was considered in Ref. '^^'^^ , while in the present paper the scattering 
on magnetic impurities in DN is taken into account. Similar to KeiM^^, we assume that the interface between the 
DN conductor and the S electrode a,t x — L has a resistance Rb while the DN/N interface at a; = has zero resistance 
and we apply the generalized boundary conditions of Ref.^^ to treat the interface between DN and S. 

We model the insulating barrier between DN and S by the delta function U{x) — H5{x — L), which provides the 
transparency of the junction Tm = 4cos^ (/ii/(4cos^ (p + Z^), where Z = 2H/vf is a dimensionless constant, </> is the 
injection angle measured from the interface normal to the junction and vp is Fermi velocity. The interface resistance 
Rb is given by 

9 
Rb = ^0 



/-7r/2 '^'?^^™ *^0^ '^ 

where Rq is Sharvin resistance R^ = e^fc|nS'c/47r^, kp is the Fermi wave-vector and Sc is the constriction area (see 
Fig. nj. Note that the area Sc is in general not equal to the cross-section area Sd of the normal conductor, therefore 
Sc/ Sd is independent parameter of our theory. This allows to vary Rd/ Rb independently of T^. In real physical 
situation, the assumption Sc < Sd means that only a part of the actual flat DN/S interface (having area Sc) is 
conducting, no matter is it a single conducting region or a series of such regions. These conducting regions are not 
constrictions in a standard sense - we don't assume the narrowing of the total cross-section, but rather that only the 
part of the cross-section is conducting. 

We apply the quasiclassical Keldysh formalism in the following calculation of the conductance. The definitions of 4 
X 4 Green's functions in DN and S, Gi{x) and G2{x), and other notations can be found in YleiM^^. The new feature 
in the present model is the spin-scattering term in the static Usadel equation-^ for Gi{x) in DN 

D^[Gi{x)^^^]+i[H-it,p„,,G^{x)]=Q, (1) 

where D is the diffusion constant in DN, H is given by 



H 



Ho 
Ho 



with Ho = eT3, and 



'^spin — -^T3Gi{x)t3 



is the self-energy for magnetic impurity scattering with the scattering rate 7 in DN. Note that magnetic impurities 
take random alignments and we average them in all directions, thus Gi{x) in our calculation is a unit matrix in the 
spin space. The Nazarov's generalized boundary condition for Gi{x) at the DN/S interface has the same form as the 
one without magnetic impurity scattering (see Ref>22i^). 

In the actual calculation it is convenient to use the standard 0-parametrization where 9{x) is a measure of the 
proximity effect in DN and is determined by the following equation 

D—^0{x) + 2i{e + i7Cos[6'(a;)]) sin[6l(x)] = 0, (2) 
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FIG. 1: Schematic illustration of the model 



One can see that introduction of magnetic impurity scattering 7 leads to modification of the effective coherence length 
in DN. In particular, switching on 7 makes function 0(x) exponentially decaying at zero energy, while 0(x) at 7 = 
behaves linearly in DN. It will be shown below that these modifications result in suppression of B in DN, as expected 
due to pair-breaking nature of magnetic scattering, which in turn leads to corresponding modifications of the subgap 
conductance. 

Finally, we obtain the following result for the electric current 
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Then the total resistance R for s-wave junction at zero temperature is given by 



with 



R 



Rb 



Rd 



dx 



ho — 



< Ibo > L Jq cosh Oim{x) 

TmAi + 2(2 - r„)A2 

2\{2-T^)+T^[gcoseL + fsm9L] \^' 



(3) 



(4) 



Ai = (1+ I cos9l P + I sin9L Hd 5 P + I / P +1) 



-4Imag[/g*]Imag[cos 6l sin 6*2] 



(5) 



A2 = Real{g(cos6lL + cos6'l) +f(sin6lL + sin6'L)}, 



(6) 



For a d-wave junction, the function Ii,q is given by the following expression 
Tn Co 



ho — 



2 I (2 - T„)(l + g+g_ + f+f_) + T„[cos0l(5+ + 5-) + sin0L(/+ + /_)] 



Co = T„(l+ I COS^i |2 + I sin^i n[| g+ + g_ f + | /+ + /- P + | 1 + /+/_ + .g+.9_ |2 + I f+g_ - g+f_ \' 

+2(2 - r„)Real{(l + g^gl + fir_)[{cos0L + cos02)(.g+ + .g_) + (sin^^ + s[n9l){f+ + /_)]} 



-4r„Imag(cos 6Il sin6'2)Imag[(/+ + ,f-){gl + g* 



g± = e/we^ — Aj_, f± = Aj-/i /Aj_ — e^ and A± = Acos2(0=Fa). In the above a, Oi„i{x) and 6'l denote the angle 

between the normal to the interface and the crystal axis of d-wave superconductors, the imaginary part of 9{x) and 
9{LJ) respectively. The conductance in the superconducting state (Js{gV) is simply related to R by (Ts(eV^) = l/R- 

It is important to note that in the present approach, according to the circuit theory, Rd/ Rb can be varied indepen- 
dently of Tm, i.e., independently of Z , since one can change the magnitude of the constriction area Sc independently. 
In other words, Rd/ Rb is no more proportional to Tav{L/l), where Tav is the averaged transmissivity of the barrier 
and I is the mean free path in the diffusive region. Based on this fact, we can choose Rd/ Rb and Z as independent 
parameters. 

In the following section, we will discuss the normalized conductance aT{e.V) = ag{eV)/af^{eV) where a!^{eV) is 
the conductance in the normal state without magnetic impurity given by aiq{eV) = ctn = ^/{Rd + Rb)- 



III. RESULTS 
A. Tunneling conductance for s-wave junctions 

In this section, we focus on the bias voltage dependent normalized conductance crT{eV) for various situations. Let 
us first focus on the relatively low transparent junctions with Z — 3 ioi various 7/A(Fig. [2J|. For Eth/^ — 1 and 
Rd/Rb = 1, the axieV) curves have a rounded bottom shape and the height of the bottom value is reduced with 
an increase in 7/ A. The height of the peak at eV = ±A is reduced with an increase in 7/ A (see Fig. I3a)). For 
Eth/ ^ = 1 and Rd/Rb = 10, the axieV) curves also have a rounded bottom structure which flattens with an increase 
in 7/A. Also the peak at eV = ±A is suppressed with the increase of 7/A (see Fig. El^b)). For small Thouless 
energy Exh/ ^ = 0.01 and Rd/Rb = 1, the conductance has a prominent ZBCP with the width given by Exh- As 
seen from Fig. O^c), the magnetic impurity scattering suppresses the peak height. With the increase of the resistance 
ratio Rd/Rb, the ZBCP transforms into ZBCD, as shown in Fig. |2Id). The magnitude of ZBCD decreases with the 
increase of 7/A, and the height of the peaks around eV/A ^ 0.04 is also reduced (see Fig. 121(d))- As seen from 
these figures, the characteristic energy range of 7 which modifies the magnitude of axieV), is determined by Exh, in 
agreement with the previous study based on the KL boundary condition^i. 

In the case of an intermediate barrier strength Z = l(Fig. EJ the magnitude of ax{eV) always exceeds unity. The 
resulting line shapes of ax{eV) for Exh/A = 1 are quite similar to the corresponding curves for Z = 3 (see Figs. 
Ela) andini^b)). For Exh/A = 1 and Rd/Rb — 1, the zero-bias value ctt(O) is independent of 7/A (see Fig|3Ia)), in 
contrast to the corresponding case shown in Fig. HJa). Another important difference from the case of large Z- factor 
is the absence of ZBCP for low Thouless energy. It is seen that for Exh/A = 0.01 a ZBCD occurs in both cases of 
Rd/Rb = 1 and Rd/Rb = 10. This conductance dip and the finite voltage peaks are fully suppressed with the increase 
of 7/A for Rd/Rb = 1 (see Fig. El^c)). On the other hand, for Rd/Rb = 10 only the peaks around eV/ IS. ^ 0.04 are 
suppressed while the magnitude of ct(0) does not depend on 7, similar to the case Z = ?, (see Fig. Eld)). The relevant 
scale of 7 is again given by the magnitude of Exh- 
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FIG. 2: Normalized conductance for Z - 



For fully transparent case with Z = 0(Fig. 0)), ax{eV) also always exceeds unity. The line shapes of ax^eV) with 
Exh/ A = 1 are similar to the corresponding curves for Z = 3 and Z = l(see Figs. 01a) and^b)). For Exh/A. — 1 
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FIG. 3: Normalized conductance for Z • 



and Rd/Rb = 1, the magnitude of axiO) is enhanced by 7/A in contrast to the corresponding cases shown in Figs, 
iga) and|3;a)(see FiglUa)). For Eth/A = 0.01 and Rd/Rb = 1, arieV) have a ZBCD. The magnitude of ariO) 
is enhanced by 7/A and the depth of the ZBCD decreases with the increase of 7/A (see Fig. EIc)). On the other 
hand, for i?T/i/A = 0.01 and Rd/Rb = 10, the magnitude of cr(0) does not depend on 7 while the finite bias peaks are 
suppressed similar to the cases of Z = 3 and Z — 1 (see Fig. Eld)). 
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FIG. 4: Normalized conductance for high transparent junctions with Z=0. 



In order to understand the above wide variety of line shapes and their relation to the proximity effect, we shall 
discuss the behavior of function 6l which is the measure of the proximity effect at the DN/S interface and determines 
the normalized local density of states by Recos9{x). At e = 0^ is always a real number even for non zero 7. First, 
we study the case oi Z — 3 and Eth/ ^ — 1 (Fig- 0) where the same values of 7/A and Rd/Rb are chosen as in Fig. |2] 
The real part of 9l has a step function like structure and it is always positive for e < A and negative otherwise. The 
absolute value of the real part of Ol decreases with an increase in 7/A. At the same time, the imaginary part of 9l 
has a coherent peak, the height of which is reduced with an increase in 7/A. For the case of Z = 3 and Eth/^ — 0.01 
(Fig. O where the same values of 7/A are chosen as in Fig. |51 the real part of 9l has a ZBCP with the width given 
by Exh- The imaginary part of 9l has a ZBCD for Rd/Rb = 1. Both the amplitudes of the real and imaginary part 
of 9l are reduced with the increase of 7/A only around zero energy in the interval of the order of Exh- 

Next we consider the case oi Z = with Eth/^ = 1 (Fig. 7) and Eth/^ = 0.01 (Fig. 8) where the same values of 
7/A are chosen as in Fig. 4. The line shapes of both Re(6'i) and lm(0i) are similar to those in Figs. |Sland|B| There 
is no clear qualitative difference between the energy dependencies of Real[Imag](0L) for Z = and those for Z = 3. 
For all cases, the magnitude of 9^ is reduced with the increase of 7 and then the proximity effect is suppressed by 
the magnetic impurity scattering within the energy range determined by Eth- In almost all cases, the magnitude of 
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FIG. 6: Real (upper panels) and imaginary (lower panels) part of Ol for Z = 3 and Eth/A. = 0.01. 

arieV) is reduced with the decrease of 9l- Only for high transparent case with not so large Rd/Rb, the decrease of 
the magnitude of 9^, i.e., the reduction of the proximity effect, can enhance the magnitude of axieV). 

In the following, we explain the wide variety of the line shapes of arieV). We consider Z = and EtkI ^ — 1 
case, where Bl has a weak energy dependence around zero voltage. For the fully transparent case with Tm = 1 , i.e., 
Z ^ 0, crT(O) can be given by 



ctt(O) 



1 + Rd/Rb 



1/ < /fco > +Rd/Rb 



(7) 



with 



< ho > 



1 + sin Ol ' 



(8) 



From this equation we find that the magnitude of ctt(O) gets close to unity under the strong proximity effect, i.e., 
when the magnitude of Rd/Rb is large. As shown in Figs. IHa) and[3^b), the magnitude of 6*^ at e = is lowered 
with an increase in 7/ A for Rd/Rb = 1. Then , according to Eqs. (7) and (8), the resulting axieV) around eV '--^ is 
slightly enhanced as shown in Fig. Qfa). For Rd/Rb = 10, the magnitude of Rd/Rb is much larger than the magnitude 
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of 1/ < ho >■ Then the 7 dependence of axiO) becomes negligible as shown in Fig. ^b). In order to understand the 
case of Z = and the small magnitude of Eth/ ^, we decompose R into Ri and R2 following the previous work™, 
where i?i and R2 are defined by 



Ri 



dx 



L Jo cosh 6'i„i(a 



and 



R2 



Rh 



Rd < IbQ > 



Fig. 9 shows that Ri has a minimum at a finite voltage which can result in a ZBCD and that R2 has a maximum 
for high transparent junctions. For a large magnitude of Rd/Rb, the effect of Ri is dominant, then the normalized 
conductance always has a ZBCD (see Figs. 9(c), 9(d) and^Jd)). Since R2 has a maximum at zero voltage (Fig. 9(b)), 
the resulting arieV) has a ZBCD as shown in Fig. 01c). 

Next we focus on the zero voltage resistance R/Rb as a function of Rd/Rb- For Z = 3, R/Rb has a reentrant 
behavior as a function of Rd/Rb due to the so called reflectionless tunneling effect^- (see Fig. 10(a)). With the 
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FIG. 9; Normalized resistance for Z = and -Brh/A = 0.01. 

increase of 7, this effect is smeared since the magnitude of 6l is reduced as shown in Fig. 11. In contrast, for Z = 0, 
where R/Rb increases monotonically as a function of Rd/Rb, the 7 dependence of R/Rb is very weak (see Fig. 10(b)). 
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FIG. 11: Real part of Ql at zero energy as a function of Rd/Rb 



B. Tunneling conductance for d-w^ave junctions 



Below we discuss the results of calculations for the d-wave case. Fig. E| shows the normalized conductance for 
Z = 10, Rd/Rb = 1, Exh/^ = 0.01, and a/ir = where a denotes the the misorientation angle between the normal 



to the interface and the crystal axis of d-wave superconductors. In this case, MARS are not formed at the interface 
of the d-wave superconductor. The origin of the ZBCP is due to the proximity effect in the DN region and the height 
of the ZBCP is suppressed with increasing 7 similar to the case of the s-wave junctions. 

With the increase of the magnitude of a the MARS are formed at the interface. The MARS contribute to the 
charge transport across the junction and leads to the formation of the ZBCP. As is seen in Fig. 13, the ZBCP 
does not depend on 7 for Z = 10, Rd/Rb = 1, Eth/^ = 0.01 and a/ir = 0.125. The similar result is obtained for 
different angle a/ir — 0.25. The reason is that MARS reduce the proximity effect in DN, therefore the influence of 
magnetic impurity scattering on the aT becomes less important. In the extreme case, a = 0.257r, the proximity effect 
is completely absent by the symmetry of the pair potential and ctt is completely independent of 7. 
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FIG. 12: Normalized conductance in a d-wave junction for Z = 10, Rd/Rb = 1, Eth/ ^ = 0.01, and a/vr = 0. 
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IV. CONCLUSIONS 



We have performed a detailed theoretical study of the conductance of diffusive normal metal / s- and d-wave 
superconductor junctions in the presence of magnetic impurities. Below the main results obtained in this paper are 
summarized. 

1. For the s-wave junctions, the proximity effect is suppressed by the magnetic impurity scattering within the energy 
range determined by the Thouless energy in DN. In this range both the real and imaginary parts of the proximity 
effect parameter, i.e., Re(0L) and lm(0i) are reduced with the increase of the magnitude of 7 for any transparency of 
the insulating barrier. 
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2. The magnitude of the normahzed bias voltage dependent conductance axieV) in the low transparent s-wave 
junctions is suppressed by the magnetic impurity scattering. On the other hand, for high transparent s-wave junctions, 
arieV) can be enhanced by the magnetic impurity scattering. 

3. In the d-wave junctions, the zero bias conductance peak formed for low transparent barriers is suppressed by the 
magnetic impurity scattering only for a '^ 0. For other misorientation angles the conductance is not sensitive to the 
magnetic impurity scattering in a diffusive normal metal. 

In the present paper, we have discussed the case where magnetic impurities are located in DN. These results can 
be also applied to the situation when the junction is in a weak magnetic field H. If the field direction is parallel to 
the junction plane, the pair-breaking rate is given by 7 = e^w^DH^/6, where w is the transverse size of the DJNt^^. 
Assuming w = 10~^r7T,, D = 10~^m^/s, A = 10~'^ey, and H = 10~^ ^ 10~^T, we estimate the pair-breaking rate 
7/A = 10^"^ ~ 10. This range of 7 corresponds to the parameters chosen in the present paper. The suppression of 
the ZBCP and ZBCD by the magnetic field was actually observed in several experiments^i^iiiiiSiiii^. The results of 
the present paper may serve as a guide to study the charge transport in the junctions with magnetic impurities or 
under applied magnetic field. 

It is also an interesting problem to study the influence of the magnetic impurity scattering on diffusive normal 
metal / triplet superconductor junctions where anomalous proximity effect is expected^^. An application of the 
present theory to the S/N/S junctions with unconventional superconductors also requires separate investigation. 

The authors appreciate useful and fruitful discussions with Yu. Nazarov and H. Itoh. This work was supported by 
the Core Research for Evolutional Science and Technology (CREST) of the Japan Science and Technology Corporation 
(JST) and a Grant-in- Aid for the 21st Century COE "Frontiers of Computational Science" . The computational aspect 
of this work has been performed at the facilities of the Supercomputer Center, Institute for Solid State Physics, 
University of Tokyo and the Computer Center. 
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